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Abstract 

The saddle point equation described by the eigenvalues of N xN Hermitian matrices 
is analized for a finite N case and the scaling relation for the large TV is considered. 
The critical point and the critical exponents of matrix model are obtained by the 
finite iV scaling. One matrix model and two-matrix model are studied in detail. 
Small N behavior for n-Ising model on a random surface is investigated. 
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1. Introduction 



The large N limit of matrix models has attracted theoretical interests in various 
fields including random matrix model, gauge theory and phase transition. Recently, 
it becomes clearer that the matrix model is closely related to the string field theory. 
Ising model on a random surface is described by a two-matrix model, in which 
the planar surface is generated in the large N limit and the up spin and down 
spin asigned by two different matrices. In this case, an external magnetic field on 
Ising spins can be introduced and the model is exactly solvable [1]. The critical 
exponents of the specific heat, the susceptibility and the spontaneous magnetization 
are obtained exactly, and they agree with the conformal field theoretical result [2]. 
Their values are different from the well known Onsager solution of 2D Ising model 
on the regular lattice. The difference is caused by the random surface, which is 
interpreted as 2D quantum gravity. The Ising spin has c=l/2 central charge due 
to fermionic character, and consequently Ising model on a random surface describes 
c=l/2 matter field coupled to 2D quantum gravity. 

The large N limit of matrix model gives planar Feyman diagrams in the pertur- 
bation of the coupling constant [3]. This planar diagrams gives us a random surface 
of the fixed genus. On each vertex of the descretized random surface, one can put 
different n-species Ising spins, which are independent on different species and only 
interact through random statistical average [4,5]. This n-Ising model on a random 
surface is particularly interesting since the matter field has central charge c=n/2 
which can be larger than one. Unfortunately, this model has not been solved yet 
for n > 2 although it is described by 2 n matrix model. Only numerical analyses 
exist in this case [4,5]. To investigate the large N limit of the matrix model, several 
techniques have been introduced. The saddle point equation [3], Schwinger Dyson 
equation, the orthogonal polynomial method [6], the double scaling method [7] and 
the renormalization group method [8] are considered. In this paper, we pursue the 
saddle point equation method for matrix models, since it is simple and it gives clear 
interpretation for the transition. Although the large N limit should be taken for 
applying the saddle point method, it may be interesting to investigate the finite N 
case, and to understand the mechanism of the phase transition together with a finite 
N scaling behavior. We hope that finite N saddle point analysis becomes comple- 
mental one to the series expansion of [4]. It is expected that our method also work 
for n-Ising case, and a small N calculation is presented for this as a preliminary 
study. 
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2. One Matrix Model 



The one matrix model has a following Hamiltonian, 

H = -TrM 2 + — TrM 4 (1) 

where M is a N x N Hermitian matrix. Denoting the eigenvalues of this matrix by 
A which are all real, the partition function is given by 

Z = /dAIIi^Ai-A^expI—S^-^E^]- ( 2 ) 

The large N limit of this partition function is evaluated by the saddle point equation 
of A. Exponentiating the measure, and differentiating the exponent by Aj we obtain 
the saddle point equation. 

-A,-|A ? + E ^^0. (3) 

There appears a critical value of g, beyond which there is no real eigenvalue of 
A. The free energy is expanded in the power series of g, 

F = T,Ck9 k (4) 

where c& behaves for large k 

c k ^A k k- 3+ ^. (5) 

The series of F is a convergent series for the large N limit, and A is —l/g c - For 
the one matrix model of (1), g c is -1/48 [3]. The exponent 7 st is called as a string 
susceptibility and becomes -1/2 for one matrix model. 

Although the saddle point equation is justified for the large N limit we apply it 
for a finite N matrix model. We approximate the large N Riemann-Hilbert integral 
equation by a finite discrete eigenvalue equation. It is easy to find the critical value 
of g c for such finite N saddle point equation and increasing number of N, g c is 
expected to become g c of the large N limit. Our aim is to develop the method 
for obtaining the critical coupling g c or cosmological constant A of (5), instead of 
calculating the partition function which requires the integration of the eigenvalue A. 
Our method is similar to Lipatov large order calculation, but for finite N the series 
of (4) becomes asymptotic and there appear k- factorial coefficient instead of (5). 
Therefore, Lipatov method [9] is not appropriate for our problem. We will show 
that this finite N saddle point method is practically useful for obtaining the critical 
value g c and the string susceptibility 7 st . 
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For example in the case of N=2, we have only two eigenvalue Ai and A2. The 
effective Hamiltonian after exponentiating the measure part, 

H ef f = \{\\ + \\) + f (A^ + A*) - ln(Ai - A 2 ) 2 (6) 

The saddle point equation is obtained by the differentiating this Hamiltonian by 
Ai and A2. There is a solution that two eigenvalues are symmetric around zero, 
Ai = — A 2 . Thus the equations reduce to one equation, 

A + 2^A 3 -^ = (7) 

This quadratic equation gives A = (— 1±(1 + %g) l l 2 )/Ag and g c = —1/8. The phys- 
ical solution appears in + sign, which becomes finite for g — > 0. Two solutions 
degenerates at g c and beyond this, there is no real eigenvalue A. The eigenvalue 
should be real since a matrix M is Hermitian. 

The obtained critical value g c = —1/8 is far from -1/48, but increasing N, we 
see that the result approaches to the correct one smoothly. We see later that the 
root singularity near g c found in N=2 case is a correct answer for N = 00. For N=3, 
considering Ai = — A3, A 2 = solution for three eigenvalues, we obtain g c = —1/16. 
When N is even number, we take a solution Ai = — Ajv, A 2 = — Ajv-i,..., and for N 
odd, we have Ai = — Ajv,...,Aat+i/2 = 0. The critical value of g c is easily obtained 
numerically by the investigation of the largest eigenvalue Ai. The obtained critical 
value g c is shown in table 1. They are indeed approaching to -1/48. 

Since the finite N scaling is expected, we discuss the critical exponent v defined 

by 

An = A 00 + ^ (8) 

where An = —l/g c for a finite N, and Aoo = 48. The exponent v is scaling exponent 
defined by 

F~(g-g c ) 2 -^f(N»(g-g c )) (9) 
Denoting the difference between Aoo and A^ by 

D N = A 00 -A N ~ N~ v (10) 

we have a ratio Rn = D N /D N _i ~ ((N — l)/N) v ~ 1 — v/N which is plotted in 
Fig. 1. The asymptotic coefficient of 1/N is estimated as v is 0.8 which is precisely 
same as the exact value of the pure gravity case with c = 0. 

It is somehow remarkable that without calculating the free energy, we obtain 
the correct values of A and the string susceptibility through a finite N saddle point 
equation up to order N=8. The string susceptibility 7 st is related to v by 

7st + -=2, (11) 
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and we have 7^ = —1/2 since v = 4/5. 

It is also interesting to study how the largest eigenvalue Ai behaves near g c for 
finite N, and to see the eigenvalue density behaves at g c . It is known that the density 
of the eigenvalue, which obeys usually the semi-circle law, becomes singular at g c . 
The density p vanishes as [3] 

p~(A-A c )§. (12) 

We find that the histogram of the density of eigenvalues for N = 8 agrees with this 
behavior. The density of eigenvalue at the critical value g c vanishes in general as 
p ~ (A - A c ) _1/7st_1/2 . 

The scaling between (A — A c ) and (g c — g) exists [10], and it is given by 

\-\ c ~(g c -g)-™. (13) 

Since j st = —1/2 for one matrix model, this scaling relation is satified apparently 
in N=2 case (7) as a solution of a quadratic equation. For the large value of N, also 
this root singularity is preserved and we checked it up to N=8. 



3. One Matrix Model with (TrM 2 ) 2 Interaction 

It is straightfoward to apply our method to one matrix model with (TrM 2 ) 2 
interaction, 

H = ^TrM 2 + —TrM 4 + -^-(TrM 2 ) 2 . (14) 

2 N N 2 V ; K J 

where g' is an additional coupling constant. When g = 0, this model becomes 

equivalent to O(N) vector model with 7 si = 1/2 and —l/g' c — 16 [11]. Recently, it 

was found that there is a critical point at g — —3/256 and g' = —9/256, where the 

string susceptibility has a positive value j st = 1/3 [12]. It is interesting to investigate 

this model by our saddle point equation, since 7 st = 1/3 is positive and clearly (13) 

can not be applied. The calculation is very easy and we have evaluated the finite N 

saddle point equation and obtained the critical values up to order N=7. The results 

of the evaluation for the critical values of g c and g' c are given in Fig. 2. The lines of 

finite N solutions smoothly converge to n = 00 solution as expected. Therefore, our 

finite N analysis works well for this model. 

The scaling relation of (10) is investigated near the critical point —l/g c = 3/256 

and —l/g'c = 9/256. For fixed —l/g' c = 9/256, we evaluated g c from the finite 

N saddle point equation. The critical value of g c is determined such that there is 
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no real solution of the eigenvalue beyond g c . In table 2, the obtained value of g c 
is represented. Using the same ratio method Rn = -Djv/-Djv-i, we estimated the 
scaling exponent v in (10). We obtained for this critical point v — 6/5. Asumming 
(11) is still valid in this case, we obtain precisely 7 si = 1/3. 

We have also verified that at a fixed g' = —1/16, near g c = 0, the ratio method 
gives v = 4/3 which leads to 7 s j = 1/2 by (11). This point at g — and g' = —1/16 
corresponds to the vector O(N) model and it is known that this model has j st = 1/2 
[12]. 

It is unexpected result that we have a correct scaling exponent v by the finite 
N saddle point equation. Usually the scaling relation of (9) and the exponent v 
are derived by the double scaling limit [7], based upon the orthogonal polynomial 
analysis for the free energy. We have restricted our investigation only on the saddle 
point equation. 



4. Two Matrix Model 

Two matrix model represents Ising model on a random surface and it is given by 
H = ^Tr(M x 2 + M\) - aTrM 1 M 2 + ^Tr(Mf + M 2 4 ) (15) 

where a is a coupling constant related to the nearest neighbour spin interaction J 
devided by the temperature kT, (5 = J/kT„ 

a = exp(-2/3) (16) 

The first and the third term of the Hamiltonian are written by the eigenvalues 
A and £ of the matrix M 1 and M 2 . For the two matrix, the Hermitian matrix is 
diagonalized by the unitary matrix and it is possible to integrate this unitary matrix. 
Then, partition function is written only by the eigenvalues with Haar measure [8], 

~ = C f UidXia^^U^iXi - A,)(& - Q exp(-U + a^A^) (17) 

where V is the first and the third term in (15), 

fo = £E(\ ? +£) + £E(*?+£) ( 18 ) 

For N=2, after exponentiating the measure part, we have an effective Hamilto- 
nian, 

Heff = (A 2 + £ 2 ) + g{\\ + & - ln[Misinh(2aA 1 £ 1 )] (19) 
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Ai + 2gA°- — - _ u * =0 (20) 



where we use Ai = — A2, £1 = — £2- We dropped the irrelevant term which is 
vanishing in the integration of eigenvalues before the exponentiation. When we take 
the solution that these eigenvalues are symmetric, Ai = £1, we have 

1 a\i 
2\~ x ~ tanh(2aA?) 

The critical value of g c is obtained from this equation. There is no real solution 
beyond g c . We obtain for a=0, g c = —1/8 and —(1 — a 2 ) 2 / g c = 16 for a — 1. The 
result that the value of (1 — a 2 ) 2 / g c at a=l becomes twice of the value at a=0 is 
consistent with the observation in the perturbation of g [A]. In Fig. 3, the value of 
g c is plotted as a function of a. Also the large N exact value of g c is given; they 
consist of two solutions [1]: the low temperature phase and the high temperature 
phase respectively given by 

_ (1 _ a y l9c = «! (21) 

1 3 a 

-d-^- " (1 ^ +a) ' 

Taking the useful analogy that the partition function of (17) is a grand canonical 
partition function with identification of g as the exponential of chemical potential 
e^, and that the canonical patition function is A = —l/g c itself, we have the free 
energy of the Ising spin on a random surface as 

F (a) = - hi A(a) (23) 

For the two-matrix model, there is a critical value of a, where the spontaneous 
magnetization vaishes. The phase described for a > 1/4 corresponds to the disorder 
phase (22) and the phase for a < 1/4 is a low temperature ordered phase (21). 

The magnetic field B for Ising spin can be introduced by the change of the 
coupling g in Vq as 

V (B) = -Tt(M 2 + Mi) + ^TrM, 4 + ^TrM 2 4 (24) 
At zero temperature a=0, the free energy F(a, B) and the magnetization M are 

F(a,B) = \ng-B (25) 

M = -dF/dB = 1 (26) 

For a finite magnetic field B, the saddle point equation gives the asymmetric 
solution, Aj 7^ £«. The low temperature phase, i.e. the symmetry breaking phase 
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also is described by the asymmetric solution. The asymmetric A 7^ £ solution is 
obtained easily, and g c is plotted by the dotted line in Fig. 3. Expanding (17) for 
small a, we have up to order a 2 , 

-£ = dW (27) 

Although it diverges at a = 1/2 when the higher order is included, instead of 
a 2 = 3/8, up to order a 2 the behavior is similar to the exact solution of (21) which 
shows the divergence at a 2 = 3/8. 

The asymmetric solution has a larger free energy than the symmetric one, and 
thus there is no symmetry breaking for N=2. The difference appears at order a 4 . 
However for a < 1/4, the difference is very small as shown in Fig. 3, and the mag- 
netization evaluated for a finite B shows effectively that there is a phase transition. 
We notice that it is necessary to subtract N(N — 1)/N terms in the expansion of a of 
exp(a X) Aj£j). Otherwise, the saddle point solution gives a wrong answer specially 
for small a. It is related to the factor of a~ N ( N ~ 1 W 2 which is divergent for a — > 0. 
Applying the saddle point equation for the two matrix model is subtle and already 
discussed in [13]. Since the first few terms in the expansion of a have no contribution 
in the integral up to order a N< - N ~ 1 ^ 2 , we safely replace exp(aX) Aj£j) term by 

E a k (£^ k /T(k + l) (28) 

fc=7V(7V-l)/2+l 

For N=5 and N=7 with the expression of (28), the symmetric solution \ = ^ 
gives the similar curves as the exact solution of the high temperature phase of (22). 
For small a, they indicate the divergence in the large N limit same as the correct 
expression of (22). Our (28) may be not sufficient for deriving the correct saddle 
point equation, since it also has terms which do not contribute to the integral. 
Indeed for even number of N, the curve of — l/g c of the symmetric solution becomes 
flat for small a. 

Since the critical value —l/g c should be same as one matrix model for a = 0, and 
the low temperature phase has a continuous curve starting at a = 0, if we obtain 
the high temperature curve which gives divergence at a = 0, it concludes that there 
is a phase transition. 

For N=7 case, we observe that the largest eigenvalue Ai has a maximum value 
at certain value of a = 0.4, and below this value Ai starts to decrease. This point 
a seems to be a critical point. This point is greater than the point where the value 
of — l/g c becomes minimum. The situation is similar to the exact solution shown in 
Fig.3. 
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More precise small a expansion using (17) is possible. Noting that the measure 
J is written by Vandermonde determinant in (17), we have invariance under cyclic 
exchange of eigenvalues and we have antisymmetric exchage between two different 
eigenvalue. The nonvanishing term in integration of (28) at order a k takes a following 
form due to these symmetries 

A^.-Aft (29) 

where l\ + ... + In — k, and U > lj for i < j. Using these nonvanishing terms in a 
small a expansion, we obtain for N=2 as 
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s- 2 +0(a*) (30) 



9c 1 - fa 

This result coincides with the previous one (27). For N=3, we find 

Since the value —(1 — a 2 ) 2 /g c at a = 1 is twice of the value at a = [4], we apply 
Pade approximation of this quantity based upon the small a expansion. Before 
making Pade analysis, we check the validity of this approximation using the exact 
expression of (21). Up to order a 2 , we have 

(I -a 2 ) 2 ia „ 2 2N 

L = 48 1 + -a 2 ) 

9c 3 

1 + M 2 , \ 

= 487— —2- 32 

where b\ and C\ are determined with the condition that at a=l we have (l+&ia 2 ) /(1 + 
Cia 2 ) = 2. We get b\ = 1/3 and c\ = —1/3. This [1,1] Pade is crude and it gives 
50.04 for a = 1/4 while the exact value at a = 1/4 is 50.63. 

Next order, up to order a 4 , [2,1] Pade approximation becomes 

a 2 ) 2 Ao (l + 8a 2 + fa 4 ) 

22 
3 



9c (TTfjj (33) 



This [2,1] Pade gives 50.36 for a = 1/4, and the result is improved. The difference 
from the exact value is 0.5 percent. Thus we see that Pade approximation is effective. 
Moreover, the second derivative of the logarithm of (33) by a shows the maximum 
at a = 0.24, which agrees with the exact critical value at a = 1/4. This means that 
we have a method of estimation of Qj C clS db specific heat peak. 
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Using the results of (30) and (31), we have Pade approximation for N = 2 and 

iV = 3 as 



= 8 1 — r~ 



(N = 2) (34) 
16(1 + a 2 ) (N = 3) (35) 



they become 8.34 and 17.0 at a = 1/4 for N=2 and N=3, respectively. 

The next order of 0(a 4 ) is calculated for N = 2 and approximated by [1,2] Pade, 

(I -a 2 ) 2 ,2,1,, l + fl« 2 

- = 8 1 + -a 2 a 4 ~ 8 7 15 n , 36 

<7 C 1 3 15 ; i + Z a 2_i| a 4 v ; 

which becomes 8.332 at a = 1/4. If we estimate the exponent v from these small 
matrix result by the ratio method (10), we obtain R 3 = D 5 /D 2 = 0.795 at a=l/4, 
which leads to v — 0.82 according to the same estimation of Fig.l. This value is 
slightly larger than the pur gravity result v = 0.8. The exact value of two matrix 
case is known as v = 6/7 = 0.853. 

The small a expansion without use of the formula of (17) is also possible. It is 
indeed necessary to develop such method for general n-Ising case, 2 n matrix models, 
since (17) is only applied to two matrix model. For example, in the N=2 case, the 
2x2 matrix is represented by 

M2=(l (37) 



b d. 

where we take Mi = diag(Xi, A2). Then Jacobian becomes 

V^i-G) 2 -^ 2 

with two eigenvalues of M2 6 and 6- Expanding a-dependent term in the exponent, 
and integrating by \b\, we have an effective Hamiltonian, 

I d \ b \ n = = exp[-(Ai + A 2 )(6 + 6) 



>/(ei-6) 2 -4|6|2 2 



+ -( Al -A 2 )x V(6-6) 2 -4|6| 2 ] 



2 



= (Al " A,) ^ 1 " &) ' n + ^(A,-A,)»({ 1 -fe)1 (39) 

where we neglected the term of order a which vanishes after integration of A, and 
also we dropped terms which vanishes for Ai = — A 2 . We have checked that the 
result of (30) is obtained by this method up to order a 2 . 
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It is possible to expand the off-diagonal elements b and trancate at the sufficient 
order both for the measure and for the exponent. The diagonal elements c and d 
are given by solving the characteristic equation in a perturbation of |6| 2 , 

c = Hi-j^- + 0{\bt) (40) 
^ = 6 + ^- + 0(|6| 4 ) (41) 

Sl - 4"2 

Jacobian is J = <9(c, <i)/<9(£i, £ 2 ), 

14" i - 42j 

After integration by these off-diagonal elements \b\ 2 , which is bounded as |6| < 
(£i —£,2)/2 , we have an effective hamiltonian written only by the eigenvalues of two 
matrices. Taking order \b\ 2 term and up to order a 2 , we obtain very close result of 
(27). Including higher expansion of \b\ 2 , the result becomes improved. This method 
can be applied for general N x N matrix. 



5. n-Ising Model 

We have considered Ising model on a random surface as two-matrix model in the 
large N limit. The extension of this model corresponds to the increasing the number 
of species or colors of Ising spin. Only same color Ising spin can interact each other. 
Denoting the number of colors by n, this n-Ising model is represented by 2 n matrix 
model [4,5]. 2 n configulations of up and down n-spins on a vertex are represented 
by 2 n matrices. 

The Hamiltonian of this matrix model is given by multi matrix model similar 
as two matrix model, in which g is a common coupling constant of TrM 4 term and 
the coefficients of TrMjMj is obtained by the inverse matrix of Boltzmann weight 
of spin interaction. 

In the previous works [4], we have calculated numerically the cosmological con- 
stant A = —l/g c , the critical value a c , and various critical exponents including the 
string susceptibility. It is desirable, however, to develop the method to calculate 
directly the cosmological constant —l/g c . The main difficulty may be that there 
is no available formula for the integration of matrix, and no systematic method to 
rewrite the Hamiltonian only by the eigenvalues of matrices. 

We apply the method which has been explained in the previous section, (38) 
or (42). Since the number of matrices increases rapidly as 2 n , we represent only 
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n=2 and n=3 Hamiltonian here. They are described by four matrices and eight 
matrices , resp ect ively. 

H(n = 2) = ^Tr(M 2 + M 2 + M 2 + M 2 ) - aTr(MiM 2 + M 2 M 3 + M 3 M 4 + M X M A ) 
- a 2 Tr(MiM 3 + M 2 M 4 ) + - Tr(M 4 + M 4 + M 4 + M 4 ) (43) 



H(n = 3) = ]- TrM? 

z i=i 

- aTr(M 1 M 2 + M 1 M 3 + M\.\ /.-, + M 2 M 4 + M 2 M 6 + ...) 

- a 2 Tr(M 1 M 4 + M : M lt + M : M- + M 2 M 3 + M 2 M 5 + ....) 

- a 3 Tr(MiM 8 + M 2 M 7 + M 3 M 6 + M 4 M 5 ) 

+ |Tr(§M 4 ) (44) 

i=l 

For N=2 case, using the representation of (37), we have for small a, 

-^rdip +0(a4) (45) 

For n-Ising case, the value of —(1 — a 2 ) 2n /g c at a=l becomes 2™ times of a=0 
[4]. Therefore, we take [1,1] Pade approximation using this condition as 

= 8 n I 7 i I 2n L21 ( 46 ) 



[1 + (-1 + 3#rry) a 

This Pade approximation gives 8.34 for n=l and 8.69 for n=2 at a = 1/4. In the 
previous series analysis [4], the critical value of —(1 — a 2 ) 2n /g c is estimated as 54.0 
at a c = 1/4 for n=2, and 59.2 at a c = 0.23 for n=4, and 63.8 at a c = 0.21 for n=6. 

It is important to find n dependence of g c for finite N matrix as (45). We have 
only discussed N=2 case, and therefore we can not make any analysis of the exponent 
v at this stage. 



6. Discussion 

In this paper, we have discussed the finite N saddle point equation and and its 
solution. We have obtained the correct scaling exponent for one matrix model and 
a modified one matrix model through finite N scaling. The exponent v has been 
determined by how finite N critical value of g c approaches to the n = 00 critical 
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value g c (10). Our result shows that even small N, the finite size N scaling works 
very well. This is rather remarkable and it may be related to the fact that the small 
g series expansion [4] gives also correct estimation of 7 si by small orders. Although 
we did not study the multicritical behavior described by TiM 21 terms, we believe 
that our method works for such cases. 

For two matrix model, we find the saddle point method is also effective although 
we restricted our analysis on small matrices. The more detailed and higher order 
analysis are necessary for this case. We have shown that small a expansion is useful 
and Pade analysis will give quite accurate numerical value for the critical value of 
g c fora finite N. 

For n-Ising model, we have only studied very small matrices. Using small a 
expansion and Pade approximation, we have discussed n dependence of the critical 
value g c . In this case also, the approach to the N = oo critical value g c seems 
smooth, and analysis is possible. This is consistent with the series expansion by 
planar diagramms [4]. There is no tachyonic instability in this model. 

The most interesting problem of n-Ising model is to find the value of 7 st for 
large n. In the previous analysis of the perturbation of g, 7 st seems increasing for 
c = n/2 > 1. For the large value of n, we need higher order calculation in the 
perturbation method [4]. Our present finite N saddle point method is complemental 
one in this respect. We are planning to study more detail for n-Ising model by 
making combined analysis with a series expansion.. 

We have considered only d = matrix models. It is also interesting to extend 
our finite N analysis for d — 1 matrix models which becomes quantum mechanical 
problems. This problem will be discussed in other place. 

Acknowledgement 

The numerical calculation is supported by a Grant-in- Aid for Scientific Research by 
the Ministry of Education, Science and Culture. 



13 



References 



[I] V. A. Kazakov, "fields, strings and critical phenomena", Les Houches 1988, 

edited by Brezin and J. Zinn- Justin, P.369, North-Holland, Amsterdam (1990). 
D. V. Boulatov and V. A. Kazakov, Phys. Lett. B187 (1987) 379. 

[2] V. G. Knizhnik, A. M. Polyakov and A. A. Zamolodochikov, Mod. Phys. Lett. 
A3 (1988) 819. 

F. David, Mod. Phys. Lett. A3 (1988) 1651. 

J. Distler and H. Kawai, Nucl. Phys. B321 (1989) 509. 

[3] E. Brezin, C. Itzykson, G. Parisi and J. B. Zuber, Commun. Math. Phys. 59 
(1978) 35. 

[4] E. Brezin and S. Hikami, Phys. Lett. B283 (1992) 203. 
S. Hikami and E. Brezin, Phys. Lett. B295 (1992) 209. 
S. Hikami, Phys. Lett. B305 (1993) 327. 

[5] C. A. Baillie and D. Johnston, Phys. Lett. B286 (1992) 44. 

J. Ambjorn, B. Durhuus, T. Jonsson and G. Thorleifsson, Nucl. Phys. B398 
(1993) 568. 

[6] M. L. Mehta, Commun. Math. Phys. 49 (1981) 327. 

[7] E. Brezin and V. A. Kazakov, Phys. Lett. B236 (1990) 144. 

D. J. Gross and A. A. Migdal, Phys. Rev. Lett. 64 (1990) 127. 
M. R. Douglas and S. H. Shenker, Nucl. Phys. B335 (1990) 635. 

[8] E. Brezin and J. Zinn- Justin, Phys. Lett. B288 (1992) 54. 

[9] L. V. Lipatov, Zh. Eksp. Teor. Fiz. 72 (1977) 411. [ Sov. Phys. JETP 45 
(1977) 216.] 

[10] M. J. Bowick and E. Brezin, Phys. Lett. B268 (1991) 21. 

[II] S. Hikami and E. Brezin, J. Phys. A12 (1979) 759. 

[12] S. Das, A. Dhar, A. Sengupta and S. Wadia, Mod. Phys. Lett. A5 (1990) 
1041. 

[13] C. Itzykson and J. -B. Zuber, J. Math. Phys. 21 (1980) 411. 



14 



Table 1. The critical value of g obtained for one matrix model with size N. 



N 


2 


3 


4 


5 


6 


7 


8 


-l/<7c 


8 


16 


21.2 


24.8 


27.4 


29.4 


31.1 



Table 2. The critical value of g at a fixed g' = -9/256 for (TrM 2 ) 2 model. 



N 


2 


3 


4 


5 


6 


7 


-l/<7c 


11.13 


25.60 


35.92 


43.29 


48.78 


52.99 
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Figure caption 



Fig. 1. The ratio method for the scaling exponent v of one matrix model. The 
ratio R N = D N / 'D N _i is plotted against l/(N + 1) where D N = 48 + l/g c . The 
slope gives v — 4/5. 

Fig. 2. The critical lines obtained from finite N saddle point solutions for g'(TrM 2 ) 2 / 'N 2 
interaction. The lines are N=2,3,...,7 from the left respectively. The dotted line is 
the critical line for N = oo in g — g' plane. The string susceptibility j st becomes 
1/3 at g = -3/256 and g' = -9/256. 

Fig. 3. The critical value of —(1 — a 2 ) 2 / g c is shown for N=2,3,5,7. The line of 
iV = oo is exact value and it has a critical point a = 1/4. Two lines of N=oo are 
expressed by (21) and (22). The dotted lines are the low temperature solutions. 
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